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^ ; Abstract 

Lattice calculations of matrix elements involving heavy-light quark bilinears are 
of interest in calculating a variety of properties of B and D mesons, including 
decay constants and mixing parameters. A large source of uncertainty in the 
p^ ', determination of these properties has been uncertainty in the normalization of the 

lattice-regularized operators that appear in the matrix elements. Tadpole-improved 
perturbation theory, as formulated by Lepage and Mackenzie, promises to reduce 
these uncertainties below the ten per cent level at one-loop. In this paper we study 
this proposal as it applies to lattice-regularized heavy-light operators. We consider 
both the commonly used zero-distance bilinear and the distance-one point-split 
operator. A self-contained section on the application of these results is included. 
The calculation reduces the value of Jb obtained from lattice calculations using the 
heavy quark effective theory. 
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1. Introduction 

For a heavy quark such as the b quark, a variety of approaches are available for 
lattice Monte Carlo calculations, including non-relativistic QCD (NRQCD) [2], the 
heavy quark effective theory (HQET) [3], extrapolation from Wilson fermion re- 
sults [4] , and formulations which interpolate continuously between these actions [5] . 
Whatever approach is used, the relationship of the lattice operators to the operators 
coming from the continuum electroweak theory must be calculated in order to make 
use of the lattice results. While these short-distance strong interaction corrections 
are in principle perturbatively calculable, in practice, the one-loop corrections are 
sometimes so large as to be of questionable reliability [6]-[9]. It is possible that 
two-loop corrections are so large as to make existing calculations very inaccurate at 
realistically attainable values of /? = Q/qq- 

It is disturbing that these large corrections exist and they seem to point to a 
deficiency in our understanding of lattice perturbation theory. Recently Lepage and 
Mackenzie [1] have suggested a reorganization of perturbation theory, the root of 
which is a nonperturbatively determined renormalization of the basic operators in 
the lattice action. These redefinitions remedy the problem of large renormalizations 
arising from lattice tadpole graphs. 

Another problem with lattice perturbation theory is that if one uses P to 
determine the perturbative coupling, aiat, one- loop perturbative corrections in 
quantities such as the mass renormalization for Wilson fermions are consistently 
underestimated. These perturbation theory problems are due to the fact that 
aiat is a poor choice of expansion parameter. For example at an inverse lattice 
spacing of 2 GeV, Lepage and Mackenzie's tadpole-improved expansion parameter 
is ay = 0.16, which is twice as large as aiat = 0.08. The reason that aiat is so different 
from ay is that aiat receives large all orders corrections from the tadpoles present in 
lattice QCD. One is led to the new expansion parameter after one redefines the basic 
fields in the lattice action and follows the tadpole improvement program referred 
to in the previous paragraph. Lepage and Mackenzie argue that the best way to 
arrive at ay is from a non-perturbative lattice determination of a perturbatively 
calculable quantity, such as the gauge field plaquette expectation value. After the 
reorganization of perturbation theory and the determination of a suitable expansion 



parameter, Lepage and Mackenzie reanalyze several lattice quantities noted for their 
poor agreement with perturbative calculations, and they find that the results of their 
program are in good agreement with the Monte Carlo results. 

The scope of the present paper is a study of the Lepage-Mackenzie prescription 
as it applies to the calculation of quantities determined using the heavy quark 
effective action, particularly fs- In the following section, we will review the effect of 
tadpole improvement on the heavy and light quark actions. To keep the explanation 
simple, discussions of one-loop perturbative corrections are postponed to Section 3. 
Then in Section 4, we determine the coupling and the scale for which the coupling 
should be evaluated as it appears in the expressions of Section 3. In Section 5, we 
complete the lattice-to-continuum matching. 

Our main results are in Tables 4a and 4b of Section 6. This concluding section 
has been made as self-contained as possible, and readers who just want to know 
how to apply our results can turn directly to it. There we incorporate the results 
of additional continuum calculations of the next-to-leading running and matching 
from the continuum effective theory at the lattice scale to the full theory at the 
scale mb [10] which is not otherwise discussed here. More detail on the application 
of those results can be found in Ref. [11]. The results are tabulated in a form in 
which they can be easily used to obtain physical predictions from lattice results. 
Throughout the paper, it is assumed that the light quark action is the standard 
Wilson fermion action. Both staggered fermion and naive fermion results correspond 
to the case r = [7] . However the conventional normalization of the Wilson fermion 
field is singular as r -^ 0, and a different convention is used in that case. This is 
discussed in the conclusions. 

2. Tadpole Improvement of the Heavy and Light Quark Actions 

Using tadpole improvement of the Wilson action for quarks on the lattice as a guide, 
one can perform tadpole improvement of the heavy quark action, and this has been 
done by Bernard in Ref. [12]. Instead of discretizing 

S= f d^xb^ {ido+gAo)b (2.1) 



as 

S = ia^ Y.n ^^(^) (^(^) ~ Uo{n-6yb{n-6)] . (2.2) 

the analog of the Lepage and Mackenzie prescription for the Wilson action is 

Stadpole-improved=ia^J2r,^^^^H^(^^ t/o (n-6) ^6(^-6) ) , (2.3) 



Uo 

where uq is defined as 

wo = (|Trt/pz„,)i/4. (2.4) 

The combination U^{x)/uo more closely corresponds to the continuum field (1 + 
iagA^(x)), than does U^{x) itself. 

Since Monte-Carlo calculations have been performed with the former of the 
two actions in Eqns. (2.2) and (2.3), we now derive the relationship between matrix 
elements determined with the two actions. Consider the matrix element 



(J°^(no, 0) J0(0)) = ^GB{no) , (2.5) 

where 

J°(n) = 6t(^)(10)7%5?(n) . (2.6) 

The two-by-four matrix (10) is the upper two rows of the projection operator 
(1 + 7o)/2. In general our notation for heavy quark fields and heavy-light currents 
follows that of Ref. [8] . As one finds by performing the b and 6^ integrations in the 
lattice-regularized functional integral in a fixed gauge field configuration {U^{n)}, 
the heavy quark propagator is just the product 

t/o(no-l,0)^..t/o(0,0)^ (2.7) 

With the tadpole-improved action, there is an additional factor of I/uq for each 
gauge field link in the product. Thus 

[^ B\'^0)\tadpole— improved ^ ^^ l^'^j 

Uq 

The B meson decay constant fs is usually extracted from numerical simulations 
by fitting GsinQ) to 

—- exp -Cnoa (2.9) 

2mB 
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As discussed in Ref. [6], for perturbative corrections, this necessitates using the 
reduced value of the heavy quark wave function renormahzation, the constant part 
of which at one-loop is, c = 4.53. 

From Eqns. (2.8) and (2.9), we see that the tadpole improvement procedure 
has no effect on the fitted value of fs- Its only effect is the change 

C^C+i^, (2.10) 

a 

that is, a linearly divergent mass renormahzation. 

Interestingly, if one uses the other fitting procedure discussed in Ref. [6] , the 
application of the Lepage-Mackenzie prescription leads one to introduce factors of 
uq which led to a final result identical to the above. Thus the reduction of c from 
24.48 to 4.53 discussed in Ref. [6] is the same reduction caused by the reorganization 
of the tadpole contributions [12]. 

We now consider the point-split operator. A variety of operators with the same 
continuum limit as in Eq. (2.6) can be constructed. Consider the set of distance-one 
bilinears which respects the remnants of the 0(3) rotational group that is present in 
the lattice heavy quark effective theory. As discussed in Ref. [8], the only operator 
in this set that one needs to consider is 

Jp%(^) = ^ E [bHn+min)\l Oh%q{n)+bHn-min-i)il 0)7%5?(n)^ . 

(2.11) 
The index i runs only over the three spatial directions. The sum of six terms has 
been chosen to respect the remnants of the 0(3) rotational group. The tadpole 
improvement procedure for this operator is to multiply it by a single factor of I/uq 
since it contains a single gauge field link. 

So far in this section, we have seen that tadpole improvement does not affect 
the extraction of /s as it is generally done in lattice Monte Carlo calculations 
using the zero-distance bilinear, and that the effect of tadpole improvement on the 
distance-one bilinear receives a contribution of I/uq. However we must still take 
into account the effect of tadpole improvement of the light quark action, and this 
will involve some additional factors. 

As conventionally defined in lattice Monte Carlo calculations the lattice op- 
erator Jg involved in calculating fs is renormalized by a factor y/2K^Z, where 
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Kc is the critical value of the hopping parameter for the light quarks. The tree 
level value of Kc = l/(8r). What Lepage and Mackenzie are advocating is a 
reorganization of perturbation theory so that the Kc factors are included in Z and the 
renormalizing factor becomes Z/(2y/r). ^ To the extent that the non-perturbatively 
calculated value of Kc agrees with its perturbatively calculated counterpart, there 
is no difference in these two prescriptions. In the next section we will see how to 
proceed with this procedure at one-loop order. 

3. Tadpole Improvement of Heavy-Light Operators at One-Loop Order 

This section will present in parallel the analysis for the zero-distance and point-split 
lattice discretizations of J5 . To summarize the previous section, we saw that rather 



than multiplying the zero-distance operator by v^k^Z, it should be multiplied by 
Z'/(2y^), and that rather than multiplying the distance-one operator by y^^/^Zps, 
it should be multiplied by Zps/{2uo^/r). In this section we will calculate and explain 
the relationship between Z and Z (and Zps and Zps) at one- loop order. 

Suppose that calculations of Z and Zps have been carried out to one-loop order 
and that the result is 

OTT J TT^ Z ("3 1 ^ 

,4 o ^'^ ) 

Zps = l + -^[ -^gps{q) + - \rv{iiaf] ■ 

oil J TX^ Z 

Suppose further that a one-loop calculation of 8kc has been performed, and that 

the result is 

1 A ^s f d Q 



1-— r-fh{q) (3.2) 

oTKc one — loop '^'^ -^ ^ 

(we will reserve the unadorned symbols t^c and uq for the non-perturbative lattice 
Monte Carlo values, and qualify them with the subscript "one-loop" when referring 
to their perturbative values). Finally, suppose that a one-loop calculation of uq 

gives, 

UO one-loop = 1 + — -^j{q) ■ (3.3) 



^ This is actually a slight generalization of the Lepage-Mackenzie prescription to 
the case r ^ 1. 



Then any quantity obtained from a lattice calculation, for which short-distance 
corrections have been correctly included at one-loop order can be multiplied by 

^'i^c one— loop \ j '-^0 one — loop \ / r, t\ 

where p and s are any small numbers, and the result remains correct at one-loop 
order. 

Without some kind of additional consideration there is no reason to prefer one 
of these one-loop results over another. In terms of the various quantities defined 
above, the tadpole improvement requires that for the Z factor of the distance-zero 
operator, p = 1/2 and s = 0. In other words the relationship between Z and Z is 



Z — ySrKcone-loopZ . (3-5) 

Thus if an expression for Z analogous to the expression for Z in Eqn. (3.1) were 
defined, we would have 

m^9iQ)+hiq)/2. (3.6) 

Similarly, for the point-split operator, p = 1/2 and s = 1, implying 



^ps '^0 one — loopy oTKq one— loop^ps \'^'' ) 

and 

9ps{q)=9ps{q)+hiq)/2+j{q) . (3.8) 

For clarity we emphasize that the relationship between Z and Z is not Z = y/8rK^Z 



and Zps = UQ^J^rKcZys- In that case, tadpole improvement would be without 
content, since we would be renormalizing by the exact same factor in both cases. 

The functions g{q), h{q), Qpsiq), and j{q) have already been calculated, and 
we now assemble the results. For various frequently appearing functions of q, we 
will use the notation Aj as defined in Ref. [13]. (Note that Ag defined in Ref. [9] 
conflicts with Aq defined in Ref. [13].) We also define A^ to be the same as Aj 



but with qo = 0. From Refs. [6] and [7], we have 

1 0(1 -q'^) l-4r2 r 

9iQ) = -Tin^ + „4 + -TTn^ - — T^ + 



4A1A2 q4 16^2 4A^^) 167r^ 
1 ^ 2 2e{l-q^) l_ 



_ 1 _^ ^ ^ _ /^U-rj _ _^. (3.9) 



_ 1 J_ .J_. e{l-q^) _ A4 + A5 , r2(2 - Ai) 



2'327r2 8Ai g4 16AfA2 4A2 ' ' 

The second and third groups are — c/2 (c is the reduced heavy quark wave function 

renormahzation discussed in Section 2) and — (/ — F)/2 (f — F is the constant part 

of the wave function renormahzation for a Wilson fermion), respectively. 

The critical hopping parameter Kc is related to the critical bare mass by the 

expression 

1 rrica 

1 + —^ • (3.10) 



8rKc 4:1^ 

Thus for massless Wilson quarks only the linearly divergent part of the correction 
to rric survives in the continuum limit. At one-loop, in terms of h{q) defined in 
(3.2), the result is [13] [14], 

^(«* = 3(2^7 -iATAjl'^'^^^^*- ''■"' 

The first term comes from the tadpole graph and has the larger value after 
integration. 

From Ref. [8], we find 

, , 4Ai - A? + 2A4 + 12r2A2 rA^f^ 
9A<1) = M + 24A,A2 ' - ^ • ^'-''^ 

where g{q) is defined in Eq. (3.9). Finally, 

i(?) = -^- (3.13) 

The integral of j(q') is determined by the coefficient of the first term in Eq. (4.1). 

The constant part of the one-loop corrections for any value of the Wilson mass 
term coefficient r can readily be obtained by numerically integrating the preceding 
expressions for g{q), h{q), QpsiQ), and j{q)- However, the scale of as{q) still has to 
be set. We turn to this in the next section. 



/? 


a-i(GeV) 


i^T^Upla,) 


ai/(3.41/a) 


Aya 


5.7 


1.15(8) 


0.549 


0.1830 


0.294 


5.9 


1.78(9) 


0.582 


0.1595 


0.198 


6.1 


2.43(15) 


0.605 


0.1450 


0.144 



Table 1. a~^, Q;\/(3.41/a) and Ay obtained from Monte 

Carlo calculations of the charmonium spectrum [16] and 

the plaquette expectation value [1] at various /3. 

4. Coupling Constant and Scale Determination 

In this section we continue with the application of the Lepage-Mackenzie pre- 
scription to determine the A-value of the coupling and the scale q* at which it 
is evaluated. Although it is in principle a higher order correction in as, a large 
source of error in the renormalization of the matrix element determining fs is 
which value of as to use. The Lepage-Mackenzie prescription for fixing the value 
of the coupling is to extract it from a non-perturbative calculation of the 1x1 
Wilson loop (i.e., the expectation value of the plaquette, Upiaq). Once it is known 
at some scale (alternatively, once its A value is known), it can be run to any other 
scale. The formula which relates the Lepage-Mackenzie perturbative coupling to the 
non-perturbatively determined (lattice Monte Carlo) plaquette expectation value is, 

- ln{^TrUpiaq) = ^av(3.41/a)[l - av(3.41/a)(1.19 + 0.025?i/) + 0{al)]. (4.1) 



The coefficient of Uf in Eq. (4.1) is for dynamical Wilson fermions with r = 1. The 
coefficient is r-dependent and is also different for staggered fermions [15]. However, 
in the application to the quenched approximation, we set n/ = 0, and the value of 
the coefficient does not enter our results. The A value for this coupling is denoted 
Ay. The results of this determination are summarized in Table 1. The values of 
a~^ are those given by the charmonium scale [16], and the values of the plaquette 
expectation value are taken from Ref. [1]. 

It remains to fix the scale g*, at which aviq) is evaluated, Lepage and 
Mackenzie propose to do that by calculating the expectation value of In(g^) in 
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the one-loop perturbative lattice correction. In equations, 

Y^Jd^-giq) 
(Hqaf) = IdU9{q)HQ^a^) (4.2) 

q*a = exp[(ln(ga)^)/2] 

Similar expressions for Yps and Yps{ln (qa)^) are defined with Qpsiq) replacing g{q). 
The only ambiguity in finding q* in either case is how to deal with the constants 
which do not arise from four-dimensional lattice integrals, and some momentum- 
space integrations involving the heavy quark field which were three-dimensional. 
Fortunately these terms are generally small relative to the lattice contributions. 
We write all such constants C as 

C = /A(^, (4.3) 

and similarly take the go dependence of the three-dimensional integrals to be 
constant. The values of the quantities just defined are tabulated in Table 2a and 2b 
for various r. These quantities were evaluated using the Monte Carlo integration 
routine VEGAS [17]. Errors on Y and y(ln(g^a^)) are at most 0{1) in the last 
decimal place. 

5. Application of Results to Lattice-Continuum Matching 

In this section, we illustrate the application of the results of the foregoing analysis 
to the lattice-to-continuum part of the matching for a couple of cases of interest. 

Now that we have Ay and a\/(3.41a~^) we use the two-loop formula with zero 
quark flavors (n/ = 0), to obtain av{q*), 

av{q) - [Po HqV^v) + (/3i//9o) Inln(gVAy))]"' (5-1) 

where (3o = ll/(47r) and (3i = 102/(47r)^. Setting n/ = is appropriate for lattice 
calculations done in the quenched approximation. The explicit dependence on the 
value of a~^ drops out of the ratio q* /Ay- Hence, the only way we have used 
results from lattice Monte Carlo calculations so far is for the expectation value 
of the plaquette. In Eq. (3.1), which gives the lattice-to-continuum matching, 
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r 


Y 


Y{\n{q'a')) 


q*a 


1.00 


-13.93 


-21.76 


2.18 


0.75 


-12.76 


-21.12 


2.29 


0.50 


-10.74 


-18.70 


2.39 


0.25 


-6.58 


-11.02 


2.31 


0.00 


2.21 


9.65 


8.88 



Table 2a. Y , (In(g^)), and q* for the zero distance 
lattice representation of J^ defined in Eq. (2.6). 



r 


l^ps 


^ps(ln(Q^a^)) 


q*a 


1.00 


-7.20 


-10.87 


2.13 


0.75 


-5.20 


-8.10 


2.18 


0.50 


-2.13 


-2.81 


1.93 


0.25 


2.98 


7.70 


3.64 


0.00 


8.65 


20.44 


3.26 



Table 2b. Y^^, (In(g^)), and q* for the distance-one 
lattice representation of J^ defined in Eq. (2.11). 

we put n = q*, thus In^ua becomes \nq*a. Thus the Monte Carlo calculation of 
the plaquette, and the perturbative calculations are sufficient by themselves to 
determine the lattice-to-continuum matching. These results are summarized in 
Tables 3a and 3b. 

It still remains to multiply Z and Zpg by the additional factors summarized 
at the beginning of Section 3, then match the continuum heavy quark effective 
theory at the lattice scale to the full continuum theory at ml . We leave this for the 
conclusions, to which we now turn. 



6. Conclusions 

In the preceding section, we illustrated the application of our one-loop results to 
obtain the factors Z and Zps, for the local and point-split heavy-light bilinears. 
From these Z's we will obtain the factors one multiplies a lattice calculation by in 
order to obtain a physical number. 

First, we must multiply by the two-loop running from q* to ml and the 
matching that takes us from the continuum heavy quark effective theory to the 
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r 


/3 = 5.7 


/3 = 5.9 


/3 = 6.1 


1.00 


0.73 


0.77 


0.80 


0.75 


0.76 


0.80 


0.82 


0.50 


0.82 


0.84 


0.86 


0.25 


0.91 


0.92 


0.93 


0.00 


1.13 


1.11 


1.11 


Table 3a. Z for various values of r and /3. 


r 


/3 = 5.7 


/3 = 5.9 


/3 = 6.1 


1.00 


0.88 


0.90 


0.91 


0.75 


0.93 


0.94 


0.95 


0.50 


1.00 


1.00 


1.00 


0.25 


1.13 


1.11 


1.10 


0.00 


1.24 


1.21 


1.19 



Table 3b. Zpg for various values of r and /3. 



continuum full theory. The continuum results are explained in detail in Refs. [10] 
and [11]. It is in this step that the determination of the scale a~^ is finally used. 
The resulting factor, Zcont^ is Q*- and a-dependent. The result after multiplying 
by this factor is tabulated in Tables 4a and 4b. Table 4a corresponds to using the 
zero-distance representation of the axial current given in Eq. (2.6), and Table 4b 
corresponds to using the distance-one representation given in Eq. (2.11). Now 
that perturbation theory has been reorganized to include tadpole corrections to 
all orders, we expect our one-loop calculation of the renormalization factor to be 
accurate to about 7%. This estimate of the magnitude of the two-loop correction 
was obtained simply by squaring the largest one-loop correction (0.27^ = 0.07) in 
Tables 4a and 4b. 

Second, we must multiply by the factors discussed in Section 3. We multiply 
Z and Zpa by l/(2y^), and 1/{2uq^), respectively. The value of uq is obtained 
from Eq. (2.4) and Table 1. We find uq = 0.861, 0.873, and 0.882 for (3 = 5.7, 5.9, 
and 6.1, respectively. 

As an example, consider the distance-zero operator at /3 = 6.1. From Table 4a, 
we find Z x Zcont to be 0.75. We multiply by l/(2y^) = 1/2 to obtain the final 
result, 0.37. It is worthwhile to compare this with the widely used value of Za 

11 



r 


(3 = 5.7 


(3 = 5.9 


(3 = 6.1 


1.00 


0.73 


0.74 


0.75 


0.75 


0.76 


0.77 


0.77 


0.50 


0.81 


0.81 


0.81 


0.25 


0.91 


0.88 


0.87 


0.00 


1.02 


0.98 


0.96 


Table 4a. Z x Zcont for various values of r and (i. 


r 


(3 = 5.7 


(3 = 5.9 


(3 = 6.1 


1.00 


0.89 


0.87 


0.86 


0.75 


0.94 


0.91 


0.89 


0.50 


1.01 


0.97 


0.95 


0.25 


1.08 


1.03 


1.00 


0.00 


1.20 


1.13 


1.09 



Table 4b. Zps x Zcont for various values of r and (3. 



of 0.8 [6], which does not benefit from tadpole improvement. To make a head-to- 
head comparison it is necessary to multiply Za = 0.8 by \^2k^ which is 0.56 at 
(3 = 6.1 [18]. The product is 0.45. Consequently, tadpole improvement results in 
a reduction of the physical value of fs by a factor of 0.37/0.45, i.e., a reduction 
of 18%. The change is similar to what one would have obtained using non-tadpole- 
improved perturbation theory with a larger, "boosted" value of a, such as a^s 
evaluated at the scale of the inverse lattice spacing. The direction of the change 
reduces the expected magnitude of the 1/m corrections. These corrections must be 
included to fit heavy quark effective field theory calculations and calculations that 
extrapolate to the 6-quark mass using propagating fermions with various moderately 
heavy masses. For a discussion of this extrapolation, see Ref. [4]. Note that in that 
work a boosted coupling is used in the determination of Z. 

As mentioned in the introduction, both staggered and naive fermion results 
correspond to the r = case [7]. The normalization convention for the Wilson 
fermion field is such that the coefficient of i()il){n) in the action is unity. This term 
is not present when r = 0, and it is necessary to choose a different normalization. 
If one normalizes the field so that in the naive continuum limit '0^'i/' has coefficient 
one, then the factor of l/(2y^) does not need to be included. The factor of 1/uq 
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remains necessary in the point-split case. 

We have apphed the Lepage-Mackenzie tadpole improvement program to the 
case of the heavy-light quark current necessary to determine fs from lattice 
calculations. We considered both the zero-distance and distance-one point-split 
lattice representations of J5 . Our analysis included a slight generalization of the 
Lepage-Mackenzie prescription to the case r 7^ 1. The results can be applied whether 
the light quark is treated as a Wilson, staggered, or naive fermion. 
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